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Abstract 



We propose the color-octet mechanism combined with the two gluon ex- 
change model for the diffractive J/ip production in hadron collisions. In the 
leading logarithmic approximation (LLA) in QCD, we find that the diffrac- 
tive J/V> production rate is related to the off-diagonal gluon density in the 
proton and to the nonperturbative color-octet matrix element of J/ip. The 
rate is found to be very sensitive to the gluon density at very small values of x 
(down to x = 0(1CP 6 )). As a result, this process may provide a wide window 
for testing the two-gluon exchange model, and may be particularly useful in 
studying the small x physics. And it may also be a golden place to test the 
color-octet mechanism proposed by solving the ip'(J/ip) surplus problem at 
the Tevatron. 

PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 
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In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
|Ij]. The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of the Pomeron and its interaction with hadrons remain a mystery. For a long time it had 
been understood that the dynamics of the "soft Pomeron" was deeply tied to confinement. 
However, it has been realized now that much can be learnt about QCD from the wide 
variety of small- x and hard diffractive processes, which are now under study experimentally. 
Of all these processes, the diffractive heavy quarkonium production has drawn specially 
attention, because their large masses provide a natural scale to guarantee the application of 
perturbative QCD. In Refs. [HH, the diffractive J/ip and T production cross sections have 
been formulated in photoproduction processes and in DIS processes in perturbative QCD. 
In the framework of perturbative QCD the Pomeron is assumed to be represented by a pair 
of gluons in the color-singlet state. This two-gluon exchange model can successfully describe 
the experimental results from HERA ||]. An important feature of this perturbative QCD 
model prediction is that the cross section for the diffractive J/ip production is expressed in 
terms of the square of the gluon density. 

So far the previous studies are focused on the diffractive processes at ep collider (pho- 
toproduction and DIS processes), we should expect that the two-gluon exchange model can 
also be used to describe the diffractive processes at hadron colliders. In this paper, we 
extend the idea of perturbative QCD description of diffractive processes from ep colliders 
to hadron colliders. Unlike in the case of the photoproduction process, at hadron colliders 
(shown in Fig.l), in the diffractive process the final state cc pair, which carry the quantum 
number of the incident gluon, is in a color-octet configuration. This color-octet cc pair must 
evolve into the physical color-singlet bound state of J/ip for the experimental observation. 
Here we use the color-octet mechanism based on the NRQCD factorization formalism || 
to describe this evolving process. In the last few years, the application of the color-octet 
mechanism has gained some important success. For example, by including this color-octet 
mechanism one might explain the ip' (J /if) surplus measured by the CDF Collaboration at 
the Tevatron ||. 

The diffractive J/ip production process in Fig.l can be realized via the color-octet 3 Si 
channel, in which the cc pair in a configuration of 3 S[ 8 ^ is produced in hard process as the 
incident gluon interacts with the proton by t channel color-singlet exchange (the two-gluon 
ladder parametrized Pomeron), and then evolve into the physical state J/ip through the long 
distance process (nonperturbative). In this long distance process some soft gluons will be 
emitted to change the color of the cc pair to coincide with the quantum number of the final 
state J/ip. The soft gluons carry little momentum and will cause little change to the final 
state J/ip spectrum in the diffractive processes. In the old color-singlet model the emitted 
gluons are hard, which will lead to a quite different momentum distribution in experiment. 
In NRQCD, the long distance evolving process is described by the nonperturbative matrix 
elements of four-fermion operators. In the process we considered in this paper, the associ- 
ated matrix element is (O 1 // ( 3 Si)^, which is also used to describe the large p t prompt J/ip 
production at the Tevatron and has been determined by several authors ||. 

Based on the NRQCD factorization formalism, as a general factorization ansatz, the rate 
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for the diffractive process gp — *> ipp^ can be factorized into the short distance part and the 
long distance part as the following form, 

\A(g P -> J/m 2 = \A(gp - (cc)[ 3 s[ 8) ]p)\ 2 x |.4((cc)[ 3 sJ 8) ] -> J/^)| 2 . (1) 

Here *A(gp — ► (cc) P^i^jp) is the short distance amplitude which can be calculated in pertur- 
bative QCD. The amplitude *4.((cc)[ 3 £'®] — > J/ip) describes the long distance evolving pro- 
cess which is related to the nonperturbative color-octet matrix element as |.4((cc)[ 3 S' 1 8 ' ) ] — > 
J/n 2 = {0iCS l ))/2Am c ^. 

For the diffractive subprocesses, gp — > J/if)p, the leading contribution comes from the 
diagrams shown in Fig.2. Due to the positive signature of these diagrams (color-singlet 
state), we know that the real part of the amplitude cancels out in the leading logarithmic 
approximation. The first two diagrams are similar to those calculated in photoproduction 
process, and the rest diagrams are new due to the existence of the gluon-gluon interaction 
vertex. From the following calculations, we can see that the diagrams other than the first 
two are needed to guarantee the gauge invariance. 

For convenience, we perform our calculations in terms of the Sudakov variables. That 
is, for the involved particles the four momenta are decomposed as, ki = otiq + fop + kiT, 
where q and p are the momenta of the incident gluon and the proton, q 2 = 0, p 2 = 0, and 
2p ■ q = W 2 = s. Here s is the total cm. energy of the gluon-proton system, a% and $ 
are the momentum fraction of q and p carried by the momentum k. kiT is the transverse 
momentum, which satisfies k^ ■ q = 0, k^ -p = 0. In our calculations, we set the momentum 
transfer t equal to zero, i.e., t = (k — k') 2 = 0. All of these Sudakov variables are fixed by 
considering the mass shell condition on the crossed lines shown in Fig.2. Take Fig. 2(a) as 
an example, the Sudakov variables are, 

k 2 M 2 + 2k 2 2k\ 

Oik — Oik' — , Pk — , Pk' — , 

s s s 

where is the mass of J/ip which is equal to 2m c . 

Following Ref. 0, the calculations are straightforward. When we perform the integral 
over the loop momentum k, the main large logarithmic contribution comes from the region 
-gr <C k^ <?C M| (i?jy is the nucleon radius) |2|. So, we calculate the amplitude as an 
expansion of k\. 

For Fig. 2(a), the imaginary part of the short distance amplitude A(g a p — > (cc)f,[ 3 Si ]p) 
is, to leading order contribution, 



Exactly, the final state of this diffractive process should include some soft gluons. However, 
becuse the soft gluons only carry little momenta, they do not change the final state kinematic 
distribution much. In this paper, we neglect the soft gluon effects in the long distance evolving 
process. 
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where F = ^-g^m c s, a and b are the color indexes of the incident gluon and the cc pair in 
color-octet 3 Si state. Factor | is the color factor. The function G(k) specifies the probability 
of finding the gluon in the proton. In the simplest three valence quark model, G = x 3. 
For Fig. 2(b), the result is, 



where the color factor is — Unlike in the case of the diffractive photoproduction processes, 
the color factors of these two diagrams (Fig. 2(a) and Fig.2(b)) are not the same. The leading 
part of the contributions from these two diagrams (which is proportional to -pr) can not 
cancel out each other. After integrating the loop momentum k, for small kj- this will lead 
to a linear singularity, not a logarithmic singularity (proper in QCD) as that in diffractive 
photoproduction process f2|. So, here in the case of diffractive process at hadron colliders, 
there must be some other diagrams to cancel out the leading part of Fig. 2 (a) and Fig. 2(b) 
to obtain the correct result. This is also due to the gauge invariance requirement. As 
mentioned above, in QCD due to the nonabelian SU(3) gauge theory there are additional 
diagrams shown in Fig.2(c)-(e) as compared with that in photoproduction at ep colliders. By 
summing up all these diagrams together, we expect that in the final result the leading part 
singularity which is proportional to pr will be canceled out, and only the terms proportional 

to p- will be retained. 

The contribution from Fig. 2(c) is, 

.^Wx^/f^Cffl, (4) 

where — \ is color factor. From Eqs.(0), (|3]), (f|), we can see that the leading part singularity 
from Fig. 2 (a)- (c) are canceled out as expected. 

By the same reason, for Fig. 2(d) and Fig.2(e), the leading part of each diagram is pro- 
portional to p-. However, their sum is only proportional to p- because the leading part is 
canceled out. Their final results is, 

^^^(-lo/f^GW- (5) 

Adding all the contributions from Fig. 2 (a)- (e), we get the imaginary part of the short 
distance amplitude, 



lmA(gp - (cc)[ 3 S[%) = F x {-JLg") J ^f w J(x',x"; k 2 T ). (6) 



Here, we rewrite the function G(k) as G(k) = f(x',x"; k\ 



i^y^T) = dG %[^f T! (7) 



where the function G(x', x"; fcf,) is the so-called off-diagonal gluon distribution function ||. 
Here, x' and x" are the momentum fraction of the proton carried by the two gluons. It is 



4 



expected that for small x, there is no big difference between the off-diagonal and the usual 
diagonal gluon densities ||. So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x', x"; Q 2 ) xg(x, Q 2 ), where x = M 2 /s. 

Finally, in the leading logarithmic approximation (LLA), we obtain the imaginary part 
of the short distance amplitude, 

lmA(gp -> (cc)[ 3 S[ 8) }p) = (-^S ab )^xg(x, Q 2 ), (8) 

where we set Q 2 = Ml. The factorization scale in the gluon density is very important 
because we know that parton distributions at small x change rapidly with this scale JT0| . In 



the literature there are various choices [p||3|,|TT|| . Here, we choose the scale to be = 2m c 
which is typically used in NRQCD factorization approach |J. 

Using Eqs.([I|) and (H), we get the cross section for the diffractive subprocess gp — ► J/ifjp, 



dajgp -> Jjjjp) \A\ 2 m^m^siQ^jOteS,) 
It |t=0 = W = "32^27 Ml )] • (9) 

We have also calculated the color-octet 1 Sq and 3 Pj subprocesses. As expected they 
do not contribute to the diffractive J/tp production at hadron colliders. That is, only the 
color-octet 3 Si contributes in this process. So the diffractive J/ip production considered in 
this paper is sensitive to the matrix element (0% ( 3 Si)"j, which is very important to describe 
the prompt J /if) production at the Tevatron ||. Therefore, the diffractive J /if) production 
at hadron colliders would provide a golden place to test the color-octet mechanism in the 
heavy quarkonium production. 

Provided the partonic cross section ([J) above, we can get the cross section of diffractive 
J /if) production at hadron level. However, as we know, there exsit nonfactorization effects 
in the hard diffractive processes at haron collisions JT^HTSf. As argued by D.E. Soper, 



these effects may be accounted by a suppression factor JF 5 . At the Tevatron, the value 
of Ts may be as small as Ts ~ 0.1 [|TJ],|r5|]. That is to say, the total cross section of the 
diffractive processes at the Tevatron may be reduced down by an order of magnitude due to 
the nonfactorization effects. So, the cross section for the diffractive process pp(p) — > J/ipp 
can be formulated as 



da(pp(p) - J/#) l69^m c a s (Q 2 ) 3 {0t( 3 S 1// 

\t=o = ^s^ — tz=- -/ dxig(x!,Q )[xg{x } Q 



dt |t_u 15 32 x 27 M§ 



■''In 



(10) 



where x% is the longitude momentum fraction of the proton (or antiproton) carried by the 
incident gluon. So, the cm. energy of the gluon-proton system is s = X\S, where S is 
the total cm. energy of the proton and proton (antiproton) system (e.g., S = (ISOOGeV") 2 
at the Tevatron). Then, x = M 2 /s = M 2 /x\S. We can see that at the Tevatron, x may 
lower down to the 10 -6 level. The lower bound of the integral variable X\ is set to satisfy 
the relation M 2 s which guarantee the validity of our calculations for this high energy 
process. 
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For the numerical calculations, we choose the input parameters as, My = 3.10GeV, 
a s (2m c ) = 0.27, (of ( 3 Si)^ = 0.0106GeV 3 . The value of the color-octet matrix element is 
taken from For the parton distribution function, we use the GRV NLO set [17]. The 
region for x is limited in GRV set to 10~ 5 < x < 1, but in our case, the value of x can reach 
the level down to 10 -6 , so we extrapolate the GRV set to these values according to their 
fitting parameters. 

In Fig. 3, we plot the differential cross section dajdtit = 0) as a function of a;i m i n for 
pp collider at the Fermilab Tevatron. (The value of the suppresion factor of J-'s is taken as 
0.1). From this figure, we see that the main contribution to the total cross section comes 
from the region of X\ > 10 -3 , which contributes 87% of the total cross section. x\ > 10~ 3 
corresponds to the region of x < 2.9 x 10~ 3 . For x\ > 10~ 2 , the contribution is over 58%, 
which corresponds to the region of x < 2.9 x 10~ 4 . For x% > 10 _1 , the contribution is above 
16%, which corresponds to x < 2.9 x 10^ 5 . We can see that the gluon density at small x 
is very important to the diffractive J /if) production at hadron colliders. Therefore, we may 
measure the gluon density at small x by observe the diffractive J /if) production. At the 
Tevatron, the gluon density can be probed down to the level of 10~ 6 through this process. 
At the LHC, the situation is even more optimistic. 

After integrating over X\, we get the total cross section for the diffractive production 
J /if) at the Tevatron, a{pp — > J/tf>p) = 66nb. Here, we assume the J /if) diffractive slope b 
to take the same value as in the photoproduction process, i.e., b = 4.5GeV~ 2 M. We have 
also estimated the gluon kx effects following ||, which will increase the total cross section 
by a factor within 5%. Our result shows that the rate of this process is more than three 
order of magnitude larger than that of large pt J / if> production in diffractive process by soft 
Pomeron exchange [IS], where the cross section at large px (> 8GeV) is of order of O.Olnb. 

However, the above results must be viewed as a rough estimate because in our calculations 
there are large theoretical uncertainties. The parton distribution function (PDF) is unknown 
to such low x region, and different sets of PDF will result in different result for the cross 
section. And the main uncertainty may come from the charm quark mass because the 
cross section sales like the eighth power of m c . Even a modest change in the charm mass 
m c will result in huge changes in the overall normalization of the cross section. Another 
uncertainty may come from the color-octet matrix element (0% ( 3 Si)^. In our estimate, we 
follow the determination of Ref. |16||. So far our formula is only to the LLA level, higher 



order corrections are neglected. Nevertheless, as a first order approximation, our result may 
provide much sense about the diffractive J/if> production rate at hadron colliders. 

Another important issue is about the nonfactorization effects of this process. In the 
above, we use the general suppression factor JF S to describe these effects. Further investi- 
gation about the nonfactorization effects is needed but is beyond the scope of this work. 
Nevertheness, we believe that the nonfactorization effects may not destroy the main results 
obtained in our paper, and the process in this paper may also be viewed as a test of the 
factorization of the hard diffractive processes in hadron collisions. 

To sum up, in this paper, we give the formula for the diffractive J / if) production at hadron 
colliders in LLA QCD by using the two-gluon exchange model. We introduce the color-octet 
mechanism to realize the color-octet cc pair evolving into J /if) meson. Though there may 
exist large theoretical uncertainties in our calculations, our results show the importance of 
diffractive J/ if) production at hadron colliders to the study of small x physics, the property of 
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diffractive process, the nature of the Pomeron and even for the test of color-octet production 
mechanism of heavy quarkonium. We hope the experimental measurement will be carried 
on in the near future. 

This work was supported in part by the National Natural Science Foundation of China, 
the State Education Commission of China, and the State Commission of Science and Tech- 
nology of China. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive J/ip production at hadron colliders in perturbative 
QCD. The black box represents the long distance process for color-octet cc pair in 3 Si state 
evolving into J /if). 

Fig. 2. The lowest order perturbative QCD diagrams for J/ip production at hadron colliders. 
Fig. 3. The differential cross section da/dt (t = 0) at the Fermilab Tevatron as a function of 
the lower bound of Xi in the integral of Eq. (|T0D . 
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FIGURES 
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Fig.l 
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